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ABSTRACT 

Heating, Ventilation and Air Conditioning (HVAC) systems con¬ 
sume almost half of the total energy use of commercial buildings. 
To optimize HVAC energy usage, it is important to understand the 
energy consumption of individual HVAC components at fine gran¬ 
ularities. However, buildings typically only have aggregate build¬ 
ing level power and thermal meters. We present HVACMeter, a 
system which leverages existing sensors in commercial HVAC sys¬ 
tems to estimate the energy consumed by individual components 
of the HVAC system, as well by each thermal zone in buildings. 
HVACMeter can be generalized to any HVAC system as it uses the 
basic understanding of HVAC operation, heat transfer equations, 
and historical sensor data to estimate energy. We deploy HVACMe¬ 
ter to three buildings on our campus, to identify the set of sensors 
that are important for accurately disaggregating energy use at the 
level of each Air Handler Unit and each thermal zone within these 
buildings. HVACMeter power estimations have on an average 44.5 
% less RMSE than that of mean power estimates. Furthermore, we 
highlight the usefulness of HVACMeter energy estimation model 
for a building fault detection application by quantifying the amount 
of energy that can be saved by fixing particular faults. 

1. INTRODUCTION 

HVAC systems consume significant energy in commercial build¬ 
ings and are naturally a focus of energy efficiency efforts. Knowing 
where the energy is being consumed at fine temporal and spatial 
(i.e. each thermal zone, air-handler unit (AHU), etc.) granularity 
is a key step towards energy optimization. Unfortunately, typical 
HVAC installation only provides aggregate electrical and thermal 
energy meters at the whole building level. While sub-metering can 
provide finer floor or subsystem level breakdowns Q, it is costly 
and still not at a fine enough spatial granularity. 

Another option to estimate energy usage in different parts of the 
HVAC system is to use simulation software such as EnergyPlus 
or DOE-2. However these are only done during the design phase 
or during extensive retrofits given the substantial cost to build ac¬ 
curate models. EnergyPlus is being extended to connect to actual 
BMS sensor points, to provide real-time, detailed energy analysis 
of HVAC systems However, these are early efforts and it takes 
considerable effort by an expert to build a calibrated EnergyPlus 
model. As a result, most building managers do not have access to 
calibrated energy models and rely on experience, aggregate meter¬ 
ing, and sensors to estimate energy use at a coarse level. 

To address this challenge of energy estimation of HVAC compo¬ 
nents, we designed ZonePac 0 using a simplified energy model of 
VAV boxes only using heat transfer equations. ZonePAC required 


detailed inspection of building architectural diagrams to apply these 
equations, and this does not port well to HVAC systems in other 
buildings. In this paper, we present HVACMeter, a system to esti¬ 
mate both the heating and cooling energy of thermal zones as well 
as AHUs, using data driven methods to apply heat transfer equa¬ 
tions without time consuming incorporation. We use HVACMeter 
on three buildings on our campus with installed building-level en¬ 
ergy meters and show that we can estimate heating and cooling 
power with an average 44.5 % less RMSE than that of mean ther¬ 
mal power. 

Facilities managers struggle to maintain HVAC fault free and en¬ 
ergy efficient 0. On our campus, across the 55 buildings managed 
by a single BMS, it is common to have over 100,000 outstanding 
faults, and more than 10,000 fault alarms are generated per day. 
While there are numerous reasons for faults being left unaddressed, 
the inability to easily prioritize them based on potential savings is 
a key limitation 0. We leverage HVACMeter for estimation of 
energy wastage due to faults detected and show that it can be in¬ 
strumental to help prioritize faults for facilities managers. Across 
three buildings we detect 19 faults with an estimated thermal en¬ 
ergy wastage of 121 MMBTU/year. 

2. BACKGROUND 

The UC San Diego campus has a central plant that generates power 
and cold/hot water circulating to over 500 buildings using campus 
wide loops. Typical building HVAC systems contain Air Handler 
Units (AHUs), as shown in Figure 0 that use this cold (hot) wa¬ 
ter, as controlled by a chilled (hot) water valve, to cool (heat) the 
incoming supply air as needed. The fiow of air and water are main¬ 
tained by fans and pumps respectively. The temperature and air 
fiow are determined by configuration points such as supply air tem¬ 
perature setpoints and supply air static pressure setpoints. Various 
sensors measure temperature of supply air, return air, static pres¬ 
sure, etc. for closed loop control. The AHUs also mix outside air 
for air quality and economizing energy use by using ambient air 
whenever possible. For each thermal zone in a building, usually 
a large room or a few smaller rooms. Variable Air Volume (VAV) 
boxes control the environment. Often there are several hundred 
VAV boxes (Figure 0 in a medium sized building. Each thermal 
zone has a thermostat to measure temperature and allow a user to 
adjust settings. The amount of cold supply air in the zone is mod¬ 
ulated using a damper, and a hot water coil can re-heat supply air 
when needed. In addition, some buildings may have specialized 
units such as dedicated chillers, exhaust fans, etc. 

We rely on the basics of heat transfer for our energy estimation. 
According to the heat transfer equations, thermal power consumed 
is proportional to the rate of flow of air/water cooled and the dif- 
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Figure 1: HVAC Air Handler Unit (AHU) 
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Figure 3: Comparison of summation of estimated VAV cooling 
power with estimated AHU cooling power 
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Figure 2: HVAC Variable Air Volume (VAV) Box 
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Figure 4: Comparison of summation of AHU cooling power 
with measured building cooling power 


ference between the initial and the final temperature of air/water. 
Hence, our goal is to estimate these physical properties from the 
sensors present for HVAC operation. However, some of these sen¬ 
sors do not directly correspond to the required rate of fiow or tem¬ 
perature measurements, and our prior work ZonePAC Q relied on 
examining building architectural diagrams to determine the rela¬ 
tionship between the sensors and physical measurements required. 
Further, the heat transfer equation cannot account for errors due 
to sensor miscalibration, losses due to inefficiency or leakage of 
air/water. HVACMeter uses historical sensor data and data driven 
methods to estimate the constants that relate sensor data to the heat 
transfer equation variables. This principle guides the estimation 
models presented in the rest of the paper. We highlight the neces¬ 
sary sensors and the minimal HVAC structural information required 
for our estimation models as well as how we can compensate for 
lack of some of this information. 

3. COOLING POWER ESTIMATION 

The cooling power consumption of a VAV box is a function of the 
supply air temperature, the return air temperature and the amount 
of air circulated. It is represented by the heat transfer equation: 

Qvav_clg — p ^ C ^ Q ^ {Tzone Tsupply) (1) 

where, Qvav_cig = VAV cooling thermal power, p = density of air 
at 20 °C, C = specific heat of air, q = supply air flow, Tzone = zone 
temperature, Tsuppiy = supply air temperature. 

Since there are no return air temperature sensors in VAVs, we ap¬ 
proximate them to the zone temperature measured by the thermo¬ 
stat. If the VAVs lack a supply air temperature sensor, they are 
approximated by the supply air temperature of the corresponding 


AHU. However, since HVAC systems can have multiple AHUs, 
knowing the mapping between a VAV and the AHU that serves it is 
needed. In some of the buildings we examined, some VAVs were 
not mapped to a thermal zone, some VAVs lacked flow sensors, and 
the AHU to VAV mapping was difficult to determine. We approxi¬ 
mated these values for these rare cases. 


We compare the total estimated power used by all of the VAVs with 
which measured for the entire HVAC system using the thermal me¬ 
ter. Note that the energy saved by the economizer is estimated sep¬ 
arately using heat transfer equations on the return air and mixed air 
temperatures. If mixed air temperature sensor is not available, it is 
estimated with outside air temperature and the economizer damper 
status. If return air temperature is not available, it is estimated as 
the average of all VAVs zone temperature. The estimated power 
also does not consider loss due to leakage or other factors. Hence, 
HVACMeter uses linear regression to estimate constant coefficients 
that better fit the total cooling power data. The final linear regres¬ 
sion formulation follows: 


vav_ _clg~\-C2^ E Qvav *{%, 


xed^TCs — Qtotal 


( 2 ) 

where, qvav = supply air flow at VAV, Qtotal = total cooling power, 
Cl, C 2 , C 3 = regression coefficients. As most AHUs do not have 
fiow sensors, we sum up the air fiow rate of each VAV that belongs 
to this AHU to obtain the total supply air flow. Figure [^compares 
the total estimated VAV cooling power (MMBTU/hour) with the 
actual value from the aggregate thermal meter. 


The estimated cooling power of an AHU is similarly a function 
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Figure 5: Comparison of summation of AHU and VAV heating 
power with building heating power measurement 


Building 

Training 

Test 

RMSE 

RMSE 

with Mean 

RMSE 

RMSE 

with Mean 


Cooling Power of AHU 

A 

0.20 

0.83 

0.26 

0.34 

B 

0.14 

0.49 

0.21 

0.52 

C 

0.042 

0.093 

0.044 

0.090 


Cooling Power of VAV 

A 

0.21 

0.83 

0.24 

0.34 

B 

0.12 

0.49 

0.18 

0.52 

C 

0.041 

0.093 

0.045 

0.090 


Heating Power 

A 

0.28 

0.50 

0.40 

0.42 

B 

0.058 

0.25 

0.074 

0.15 

C 

0.019 

0.024 

0.022 

0.02 


(Unit: MMBTU/Year) 
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Figure 6: VAV fault: Air flow does not match its setpoint 
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Figure 7: Economizer damper fault: Estimated and measured 
mixed air temperature do not match 


Table 1: Root Mean Squared Error of our estimation models 
against measured cooling and heating power as measured by 
whole building thermal meters versus error introduced if mean 
power was used as an estimate. 

of the supply air temperature, the mixed air temperature and the 
supply air flow. The summation of the cooling power of all AHUs 
in a building should match the measured total cooling power. We 
use following linear regression: 

Qahu — (Tmixed Tgupply^ * ^ ^ Qvav ( 3 ) 

C 4 * ^ ^ Qahu T — Qtotal (4) 

where, Qahu = AHU thermal power, Tmixed = mixed air temper¬ 
ature, Tsuppiy = supply air temperature, qvav = supply air flow at 
VAV, C 4 , C 5 = regression coefficients. When if Tmixed < Tguppiy, 
then AHU is in heating mode. We Alter the data for Qahu for equa¬ 
tion for AHU cooling mode only. Figure [^compares the sum of 
estimated AHU cooling power with measured cooling power. 

4. HEATING POWER ESTIMATION 

Unlike cold water, hot water is used by both VAVs and AHUs as 
shown in Figures andj^ The heating power used by an AHU 
is the same as equation ^ except that the data is now Altered for 
AHU heating mode and the sign has changed. For VAVs, however, 
there is no supply air temperature sensor available. The only points 
which indicate use of heating is “heating valve command” which 
controls the flow of water in the heating coil. The supply water 
temperature sensor is available, but the return water temperature 
is not, so we cannot apply the heat transfer equation as we do for 


AHU. We approximate the heating power at VAV as follows: 

Qvav_htg — (Tgupplywater Tgupplyair^ * Qvav * Hyalve ( 5 ) 

where, Qvav_htg = VAV heating power, Tsuppiywater = hot water 
supply temperature, Tsuppiyair = supply air temperature, Hyaive 
= heating valve command. This approximation is justifled because 
VAVs heating valve command is linearly proportional to the amount 
of hot water used for HVAC systems we examined. 

We again use linear regression as following: 

Cq * ahu + C 7 * vav_htg Cs — Qtotal ( 6 ) 

where, = regression coefficients. Figure compares 

total measured heating power with the sum of estimated heating 
power. Note that Qahu is estimated using mixed air temp, and thus 
the savings due to economizer heating are already included. 

Table compares the estimated values with the measured power 
consumption. We compare the RMSE of our estimate with the 
RMSE if just mean power were used as an estimate. The results 
show not only a good match of the estimate to actual energy con¬ 
sumption over time, but also as summarized by Table the RMSE 
of HVACMeter is less than RMSE of mean power by 53%, 56%, 
and 27% for AHU cooling power, VAV cooling power, and AHU-i-VAV 
heating power respectively. The variation in energy estimate across 
buildings is signiflcant and subject of our continuing study. This 
estimation has another important use in identifying faults that we 
discuss next. 

















































































































Index 

Rules 

Possible Faults 

Detected Faults 

(Number) 

Energy Loss 
(MMBTU/Year) 

1 

CoTTClCLtiOTli^M A, measured ^Tm A, estimated) ^ Cth,l 

Economizer damper broken 

3 

55.56 

2 

RMSPEi^Qyav •) QSvav) ^ Cth,2 

VAV damper leaking or stuck 

7 

49.40 

3 

UTlOCCUpied & (Tzone Tiipp^yiifYiif,), Qvav ^ 1.1 • Qmin 

Configuration error 

7 

14.58 

4 

When Cvaive = 0, MPE{Tsa,ahu, Tma) < Cth,3 

AHU cooling valve leaking 

2 

1.349 


T=Temperature, qs = Supply Flow Setpoint, SA=Supply Air, MA = Mixed Air, 

Cvaive=^oo\mg Vulve Command, MPE = Mean Percentage Error, C'th=Threshold q=Supply Flow, 

AHU=Air Handler Unit, VAV=Variable Air Volume, RMSPE=Root Mean Square Percentage Error, 

Table 2: Fault Detection Rules and Result 


5. FAULT PRIORITIZATION 

We detect specific types of faults in AHUs and VAVs that are known 
to facility managers, but are not detected by fault detection tech¬ 
niques used by BMS. We use rule based detection method s |[5| [7), 
as they are supported by many fault detection frameworks and 
can be designed to detect any type of fault. Table [^lists the type of 
faults we detect and a summary of the rules used for detection. 

The generic fault detection methodology we use is to estimate the 
ideal values for a parameter, and compare it with its measured 
value. For some faults, the estimation is simple, as there points in¬ 
dicating the setpoint for some sensors. For example, measured sup¬ 
ply air fiow should be close to its setpoint. For other faults, we use 
sensor data and knowledge about HVAC working for detection. For 
example, we can check if the economizer damper is working cor¬ 
rectly if the mixed air temperature changes with changes in damper 
position. We estimate the mixed air temperature by using mea¬ 
sured return air temperature, outside air temperature, and econo¬ 
mizer damper position, which determines the proportion of outside 
air and return air in the mixer. We can compare this estimated value 
to the measured mixed air temperature for detecting faults. How¬ 
ever, it is common for temperature sensors to develop minor drifts 
over time, so to detect damper faults more accurately we check the 
correlation between the measured and estimated mixed air temper¬ 
ature. When they are not correlated, either one of the sensors is 
incorrect, or the outside air damper is damaged. Figure]^ shows an 
example of a fault detected using this metric, and we can see that 
the observed mixed air temperature does not change as expected 
with changes in economizer damper position. 

To reduce false positives, and to ensure we detect faults which 
cause significant wastage or discomfort, we use conservative thresh¬ 
olds in our rules and only identify faults that persist for extended 
periods of time. We ignore transient and minor faults as facilities 
management does not have enough manpower to fix them across 
55 buildings. Figure shows an example of supply air flow fault 
detected using our rules, and we can observe that the air flow is 
almost zero even when the setpoint is set to 300 cfm. We detected 
total of 19 faults across three buildings (Table|^. 

Most importantly, to prioritize the faults that cause egregious wastage, 
we estimate the energy impact of faults using the simplified HVACMe- 
ter energy model we presented in Sections and We estimate 
energy loss for those faults whose ideal values can be estimated by 
comparing it to the energy consumption in faulty operating condi¬ 
tions. We calculate savings based on at least one week of data and 
extrapolate the savings for a year based on correlation with out¬ 
door weather conditions. Table |2] summarizes the faults detected 
and their energy wastage. For the 19 faults detected, we estimate 


energy wastage of 121 MMBTU/year using HVACMeter. 

6. DISCUSSION AND FUTURE WORK 

We have shown that it is possible to disaggregate HVAC thermal 
power consumption into individual zones and AHUs. Though our 
estimated values are not as accurate as those provided by detailed 
simulation engines like EnergyPlus, HVACMeter provides us in¬ 
sight into energy fiows within the HVAC system with minimal ef¬ 
fort. With HVACMeter, it is easy to identify zones which consume 
high amount of energy and stress the AHUs, as well as identify 
AHUs which have abnormally high energy use due to faults, and 
then optimize their operation. 

For data inference applications like HVACMeter, HVAC systems 
need to have a standardized ontology and historical sensor data, 
both of which are not common in many buildings. If they are avail¬ 
able, innovative methods can be developed to exploit existing in¬ 
formation to extract insights about HVAC operation. For exam¬ 
ple, AHUs across multiple buildings could be compared to check if 
economizers are functioning correctly. In addition, such data anal¬ 
ysis will identify which sensors are essential for data driven mon¬ 
itoring of HVAC systems, and how to compensate for faulty and 
miscalibrated sensors. 
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